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The Entropy of the Quark-Gluon Plasma 
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The entropy of the quark-gluon plasma can be calculated from QCD using (approxi¬ 
mately) self-consistent approximations. Lattice results for pure gauge theories are accu¬ 
rately reproduced down to temperatures of the order of 2.5Tc. Comparisons with other 
approaches to the thermodynamics of the quark-gluon plasma are briefly discussed. 


1. Introduction 

As we were reminded, the Bielefeld meetings which took place in 1980 and 1982 have 
played a decisive role in our held. It should be added that, since then, the Bielefeld group 
has played a leading role in the study of QCD thermodynamics with lattice calculations. 
There is one particular result of that time that I wish to recall here because it is central 
to the topic of my talk. This result is that of the lattice calculation of the energy density 
of SU(2) pure gauge theory [^. It is reproduced in Fig. 1 where one sees clearly the 
phase transition and the approach to the ideal gas limit at high temperature conhrming 
expectations based on asymptotic freedom. 

Modern results for the SU(3) equation of state ^ are displayed in Fig. 2. The talk 
will focus on the slow approach to the free gas limit of the thermodynamical functions at 
high T. The question that I wish to address is whether one can understand the deviation 
from the ideal gas behaviour with weak coupling calculations. Beyond this question is, 
as we shall see, another one related to the relevance of the notion of quasiparticles in the 
description of the quark-gluon plasma. 

The motivation for considering the quark-gluon plasma as a weakly coupled system is 
of course asymptotic freedom, which leads us to expect that the effective coupling g to 
be used in thermodynamical calculations should be small if the temperature T is high 
enough: 
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with typically /i ~ 27rT. But we know that, even in cases where the coupling is small, 
strict perturbation theory cannot be used. Technically, infrared divergences occur in high 
order calculations, and various resummations are needed to get meaningful results. 

Some of the difficulties of perturbation theory are already visible on the lowest orders 
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Figure 1. Energy density of a SU(2) gluon plasma normalized to that of an ideal gas of 
gluons, as a function of the temperature. From Ref. [1]. 


which, for SU(3) gluons, read: 


P 

'Po 




+ 30 



with Pq the ideal gas pressure: 
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The large coefficient of the term makes its contribution to the pressure comparable 
to that of the order when g ^ 0.8, or ag ^ 0.05; larger coupling makes the pressure 
larger than that of the ideal gas. Now it is worth emphasizing that the term of order g^ 
emerges from an infinite resummation (strict perturbation theory at finite order would 
lead to a polynomial in g'^)] it is the leading term of the resummed expression when g is 
small. We shall later argue that the underlying strategy (that of performing an infinite 
resummation which is then re-expanded in powers of g) may not be used to extrapolate to 
large coupling. Physically the need for resummation arises from the existence of collective 
excitations in the system, whose properties are not well captured by perturbation theory. 

A lot of efforts have been devoted in the recent years to push perturbative calculations 
of the pressure to the highest calculable order, namely order g^ [^]-|^. The contributions 
of order g'^ and g^ do not resolve the difficulties met at order g^: the values of the pressure 
obtained by adding successively these high order contributions oscillate wildly, and reveal 
a strong dependence on the renormalization scale. Attempts have been made to construct 
smooth extrapolations based on the first terms of the series, using Fade approximants 
or Borel summation techniques p,^]. The resulting expressions are indeed smooth 
functions of the coupling, better behaved than polynomial approximations truncated at 
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Figure 2. The thermodynamical functions in SU(3) gauge theories. From Ref. [2]. 


order or lower, with a weak dependence on the renormalization scale. However these 
techniques, which can be powerful, offer little physical insight. 

A more appealing physical picture is provided by quasiparticle models in which one 
assumes that the dominant effect of the interactions can be incorporated in the spectral 
properties of suitably defined quasiparticles. It is useful to review briefly the simplest 
version of such models, limiting ourselves here to bosonic quasiparticles []^. One assumes 
that the system is composed of non interacting massive quasiparticles with energies Ek'. 


El = k^ + m^{T) 


(4) 


with the mass m some function of the temperature. The entropy density is that of the 
ideal gas: 

r d?k 


s = 


(27r) 


[(1 + Nt) log(l + Nt) - Nt \ogNt ], 


where 
Nk = 




(5) 


( 6 ) 


The pressure P and the energy density e are then given by the corresponding ideal gas ex¬ 
pressions, to within a function B{T) adjusted so as to satisfy thermodynamical identities. 
That is, one sets: 


P = -T 


d^k 


log (l - - B{T) 


(7) 
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iVfcEfc + s(r). 


( 8 ) 


Such a parametrization obviously fulfills the identity e + P = Ts. The function B{T) is 
then determined by requiring that 


s = 


(W 

dT’ 


(9) 


Such a model has been shown to provide a quite good basis for a fit to lattice results 
Recent extensions of the model can be found in Ref. . 

Finding a microscopic justification for this quasiparticle picture is the main motivation 
of the work that I shall present, summarizing results obtained in collaboration with E. 

]I3|-[I^. Our analysis is based on the entropy, for which the simple 
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formulae above already suggest why it is simpler than the other thermodynamical func¬ 
tions: in both the pressure and the energy density there are interaction contributions, 
which cancel out in the entropy, leaving an expression in which all interaction effects 
are summarized in the properties of quasiparticles (in the simple model, the temperature 
dependent mass m{T)). We shall see later under which conditions this simple picture is 
valid. 

At this stage, it is useful to have a more general characterization of the effects of weak 
interactions in a quark-gluon plasma. Let me then go through a short digression about 


the basic scales and degrees of freedom of such a system when the coupling is small [17 


2. Scales and degrees of freedom. Hard thermal loops 


In the absence of interactions, the quark-gluon plasma is a gas of massless particles 
with energy E^. = k. Small interactions have little effects on most particles which have 
momentum k ^ T, but they strongly modify the propagation of low momentum modes. 
Since the coupling to gauge fields occurs typically through covariant derivatives, = 
dx+igA{x), the effect of interactions on particle motion depends indeed on the momentum 
of the excitation and the magnitude of the gauge field. A measure of the strength of 
the gauge fields is obtained from the magnitude of their thermal fluctuations, that is 
A = {A‘^{t, x)). In equilibrium {A^(t,x)) is independent of t and x and given, in the 

non interacting case, by: 


(A^) 


d^k Nk 


( 10 ) 


Here we shall use this formula also in the interacting case, assuming that the effects of the 
interactions can be accounted for simply by a change of (a more complete calculation 
is presented in Ref. 0). 

For the plasma particles Ek = k ^ T and {A‘^)t ~ The associated electric (or 
magnetic) field fluctuations {{dA)^)T ~ are a dominant contribution to the plasma 
energy density. As already mentioned, these short wavelength, or hard, gauge field fluctu¬ 
ations produce a small perturbation on the motion of a plasma particle with momentum 
/c ~ T, since gA^ ~ gT k. However, this is not so for an excitation at the momentum 
scale k ~ gT, since then the two terms in the covariant derivative and gA^ become 
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comparable. That is, an excitation with momentum gT is non perturbatively affected 
by the hard thermal ffuctuations. And indeed, the scale gT is that at which collective 
phenomena develop. The emergence of the Debye screening mass ~ gT is one of the 
simplest examples of such phenomena. 

The ffuctuations at the soft scale gT <C T can be described by classical fields. In fact 
the associated occupation numbers are large, and accordingly one can replace by 
T/£'fc ~ l/gf in eq. (0). Introducing an upper cut-off gT in the momentum integral, one 
then gets: 

{A^)gT ~ y ^ ~ (11) 

Thus AgT ~ y/gT so that gAgx ~ g^^'^T is still of higher order than the kinetic term 
dx gT. In that sense the soft modes with k ~ gT are still perturbative, i.e. their self¬ 
interactions can be ignored in a first approximation. Note however that they generate 
contributions to physical observables which are not analytic in g"^, as shown by the example 
of the order g^ contribution to the pressure (or as in eq. fpTj) above). 

At the lower momentum scale k ~ g^T, unscreened magnetic ffuctuations play a domi¬ 
nant role. The contributions of such ultrasoft ffuctuations is 

(AVr~Tyy ^d'fcl ~ g^T\ (12) 

so that gAg 2 rp ~ g‘^T is now of the same order as the ultrasoft derivative ~ g^T: the 
ffuctuations are no longer perturbative. This is the origin of the breakdown of perturbation 
theory in high temperature QCD. A more detailed analysis reveals that the ffuctuations 
at scale g‘^T come from the zero Matsubara frequency and correspond therefore to those 
of a three dimensional theory of static fields. 

Having identified these various scales, we note now that the thermodynamical functions 
at high temperature are dominated by hard degrees of freedom. The quasiparticle picture 
that we shall develop takes into account consistently the contributions of the hard degrees 
of freedom together with those of the long wavelength collective excitations. We shall 
comment later on the techniques based on dimensional reduction which allows a treatment 
of the ultrasoft contributions. 

There exists a well developed effective theory for the soft collective excitations that 
we briefly outline now. These soft excitations can be described in terms of average fields 
which obey classical equations of motion. In QED these equations are Maxwell equations: 

d,F'^^=jL, + fext (13) 

with a source term composed of an external perturbation and an extra contribution 
referred to as the induced current. The induced current, which summarizes the effects 
of the hard degrees of freedom, can be obtained using linear response theory. It is of the 
form: 

j i*yn%(x-y)A-(y), (14) 

where A'^ijj) is the total gauge field (including the induced potential) and the (retarded) 
response function — y) is also referred to as the polarization tensor. In leading 
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order in weak coupling, this polarization tensor is given by the one-loop approximation, 
and can be written as n(a;,p) = g‘^T‘^ f{uj/p,p/T) with / a dimensionless function. Fur¬ 
thermore, the relevant kinematics is such that the incoming momentum is soft while the 
loop momentum is hard. Thus, in leading order in p/T ~ gf. If is of the form g^T‘^ f {uj/p). 
This particular contribution of the one-loop polarization tensor is an example of what has 
been called a “hard thermal loop” . 

In a non Abelian theory, linear response is not sufficient: constraints due to gauge sym¬ 
metry force us to take into account specihc non linear effects. The relevant generalization 
of the Yang-Mills equation reads |T^J^ : 

1 ___ (15) 


= n J 


\ _ , ^ 

' r\ ^ fiup-^b-" -^C 


where the induced current in the right hand side is non-linear: when expanded in pow¬ 
ers of it generates an inhnite series of hard thermal loops (self-energy and vertex 
corrections). 

In the next section, we shall consider various reorganizations of the perturbative expan¬ 
sion inspired by the previous analysis and which try to accommodate the non perturbative 
features associated with the various scales. 


3. Reorganizing pertnrbative expansion 

As we have seen, one of the main effects of the thermal fluctuations is to generate a 
mass for the soft modes. This is not peculiar to gauge theories, but occurs also in the 
simpler case of a scalar held. In this case, the mass generation can be taken into account 
non perturbatively by a very simple reorganization of the perturbative expansion, leading 


to the so-called “ 
1 


C = Co 

= Cq + C^^^ , 
with C'o = Cq — (l/2)m^0^ 


screened perturbation theory” 
1 


m- One writes: 


+ Cmt 


(16) 


The mass m ~ gT which appears here is for instance that 
obtained in leading order from the hard thermal loop. A perturbative expansion in terms 
of screened propagators (that is keeping the screening mass as a parameter, i.e. not as 
a perturbative correction to be expanded out) has been shown to be quite stable with 
good convergence properties. However, m makes the propagators explicitly temperature 
dependent, and the ultraviolet divergences which occur in high order calculations be¬ 
come temperature dependent. Although we know that temperature dependent inhnities 
should eventually cancel out, the systematics of such cancellations is not immediately 
transparent. 

In the case of gauge theory, the effect of the interactions is more complicated than just 
generating a mass. But we know how to determine the dominant corrections to the self¬ 
energies. When the momenta are soft, these are given by the hard thermal loops discussed 
above. By adding these corrections to the tree level Lagrangian, and subtracting them 
from the interaction part, one generates the so-called hard thermal loop perturbation 
theory 12^^!: 

. (17) 


C — Cq c 


HTL 
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HTL 
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The resulting perturbative expansion is made complicated however by the non local na¬ 
ture of the hard thermal loop action, and by the necessity of introducing temperature 
dependent counter terms. Also, in such a scheme, one is led to use the hard thermal loop 
approximation in kinematical regimes where it is not justihed. 

Another approach that I wish to mention at this stage is that based on dimensional 
reduction. This approach puts emphasis on the non perturbative, very long wavelength 
fluctuations. It has been shown in Refs. how, in principle, an effective theory could 

be constructed to deal with this particular problem by marrying analytical techniques 
(to determine the coefficients of the effective theory) and numerical ones (to solve the 
non perturbative 3-dimensional effective theory). The resulting effective theory is a 3- 
dimensional theory of static fields, with Lagrangian: 

A// = + SC, (18) 

with Di = di — igy/TAi. This strategy has been applied recently to the calculation 
of the free energy of the quark-gluon plasma a high temperature This technique 

of dimensional reduction puts a special weight on the static sector (it singles out the 
contributions of the zero Matsubara frequency), and a major effort is devoted to the 
calculation of the coefficients of the effective Lagrangian (which contain the dominant 
contribution to the thermodynamical functions). 

In contrast, the approach that we have developed keep the full spectral information 
that one has about the plasma excitations. In a sense, it is close to the approach based on 
hard thermal loop perturbation theory. It differs in that we focus on the simplest of the 
thermodynamical function, namely the entropy, for which indeed interesting cancellations 
occur which allow us to bypass some of the difficulties of hard thermal loop perturbation 
theory. 


4. Propagator renormalisation. Skeleton expansion 


As already stated, we want to incorporate full spectral information about quasiparticles 
in thermodynamical calculations, and we are therefore led to carry out a propagator 
renormalisation. Although in gauge theory vertex renormalisation should be carried along 
with propagator renormalisation, we shall see here that special kinematical conditions 
allow us to develop useful approximations without having to do so. 

The formalism that we shall rely on was developed long ago in the context of the 
non relativistic many-body problem. It is based on an expression for the pressure as a 
functional of full propagators, first written out by Luttinger and Ward p3. Systematic 


developments using either functional Legendre transforms or diagrammatic formalisms 
were presented later by De Dominicis and Martin |^. In implementing these methods 


to field theory one meets extra difficulties related to ultraviolet divergences. I now briefly 
outline the general method, writing explicitly only formulae appropriate for a scalar field 
(similar ones hold for QCD). 

The free energy T can be written as the following functional of the full propagator D\ 


r\D] = - lTVn£> + 4|£>], 


(19) 











where <h is the sum of all the two-particle irreducible (skeletons) diagrams. The self-energy 
n is related to the propagator by Dyson’s equation 


= D„-‘ + n, 


( 20 ) 


and Dq is the bare propagator. A remarkable feature of the functional ( p!9[) is its station- 
arity property, i.e., 6J^[D]/6D = 0 when 


(5<h 

n[D] = 2 —. 

^ ^ 5D 


( 21 ) 


Together Eqs. 0 and m form a self-consistent set of equations which determine the 
full propagator. A major observation by Baym is that self-consistent approximations can 
be dehned for any selection of skeleton diagrams in $; furthermore such “$-derivable” 
approximations have the property to respect the (global) symmetries of the hamiltonian 

II- 


The stationarity property of the free energy entails important simplihcations in the 
calculation of the entropy. Indeed we have 


5 = - 




dJ^ 

W 


D 


( 22 ) 


where the last step uses the fact that the temperature dependence of the propagator can 
be ignored. Explicitly one gets: 


5 = 


d^k dN{uj] 


ImlogD ^{u!,k) + j- 


(27r)4 dT 
where N{uj) = l/(e^‘^ — 1), and 

a(T<h) 


d^k dN{uj) 
(27r)4 dT 


Imn(a;, k)ReD{uj, k) -|- S' (23) 


5' = — 


dT 


r d^k dN{uj) 

J (27r)4 dT 


Rell ImD. 


A further simplification occurs when one restricts the choice of skeletons to low order ones 
Thus, for the two-loop skeletons in QCD, S' = 0. 


5. The 2-loop entropy 

At two loop in the skeleton expansion, the entropy takes then the simple form 

H: 

^ ^ If 

The simplihcations discussed above have led to important cancellations leaving for the 
entropy an expression which is effectively a one-loop expression, thus emphasizing the 
direct relation between the entropy and the quasiparticle spectrum. Residual interactions 
start contributing at order 3-loop. This expression is also manifestly ultraviolet-hnite. 

Albeit simple, the 2-loop expression for the entropy is nevertheless difficult to use as 
it stands. Recall that the propagator and self-energy in this equation are solution of a 
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self-consistent eqnation which is in general difficnlt to renormalize and to solve. Besides, 
in the case of gange theories, the scheme does not fnlly respect gange symmetry. We have 
therefore looked for additional approximations which exploits the simplicity of the scheme 
(the one-loop featnre of the entropy formnla) and at the same time allows for practical 
calcnlations. The main idea is to hnd gange invariant approximations for the self-energy by 
exploiting kinematical approximations snch as those which lead in particnlar to the hard 
thermal loops. The entropy is then calcnlated exactly with no fnrther approximation. One 
constraint that we impose on onr approximate 11 is that, once inserted in the formnla (p4|), 
the resnlting valne of S is pertnrbatively correct np to order (which is the maximnm 
pertnrbative accnracy that one can achieve with 2-loop skeletons). 

In the regime where the loop momenta are soft, we can nse as an approximation for 11 
the corresponding hard thermal loop (HTL), that is 


P ^ gT T^soft ~ ^HTL 

For hard momenta on the other hand, the hard thermal loop is no longer valid: corrections 
to hard particle dispersion relations dne to their conpling with soft modes need to be 
taken into acconnt. This can be estimated nsing nsnal HTL pertnrbation theory. A 
simplihcation occnrs however because we need such corrections only near the light cone: 

cn, p ~ T Uhardiuj'^ ~ P^) 


In fact we shall find convenient to dehne two successive approximations. In the first 
one, called the HTL approximation, we use H = Hhtl at all momenta. This is a hne 
approximation for soft momenta, but it is also a decent approximation at hard momenta 
because the self-energy is needed only in the vicinity of the light cone where the hard 
thermal loop approximation coincides with the full one-loop result. We call Shtl fhe 
resulting expression for the entropy. This expression fully accounts for the perturbative 
contributions of order but for only 1/4 of that of order g^. In the next-to-leading 
approximation one uses Usoft = ^htl and Hhard = ^htl + ^B, where (fH contains the 
coupling to the soft modes (I should emphasize at this point that the construction of our 
next-to-leading approximation involves several subtle steps which I am skipping here but 
which are detailed by A. Rebhan in his contribution |^). We call Snla the corresponding 
entropy; it takes fully into account all perturbative contributions of order and p^. 

As an illustration of the quality of the results which we can obtain, I present in Fig. 3 the 
entropy of pure SU(3) gauge theory normalized to the ideal gas entropy. The agreement at 
large T (T ^ 2.5Tc) is quite good. Note also that in going from one level of approximation 
to the next (i.e. from Shtl to Snla), the changes are moderate, in contrast to what 
happens in ordinary perturbation theory. This reflects the stability of the present scheme. 
It also points to the fact that the contribution of the soft collective modes is indeed small: 
the fact that they give a seemingly large contribution at order p^ in perturbation theory 
is just an artifact of the truncation of a resummed expression at hnite order. 


6. CONCLUSIONS 

Our results conhrm that for temperatures larger than about 3Tc, the thermodynamical 
functions of the quark-gluon plasma, in particular its entropy, can indeed be interpreted 
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Figure 3. The entropy of pure SU(3) gauge theory normalized to the ideal gas entropy 
SsB- Full lines: Shtl- Dashed-dotted lines: Snla- For each approximation, the two lines 
correspond to the choices jl = ttT and fl = AnT of the running coupling constant as{fi) 
in the MS renormalization scheme. The dark grey band represents the lattice results of 
Ref. [2]. 


as those of weakly interacting quasiparticles. The dominant effect of the interactions is 
to modify the spectrum of these quasiparticles. As the coupling grows the quasiparticle 
properties are non perturbatively renormalized but their mutual interactions remain weak. 

Thermodynamic functions are dominated by hard degrees of freedom. While long range 
correlations may survive in the quark-gluon plasma at very large temperature, the present 
picture suggests that such correlations do not contribute much to the thermodynamics. 
We hnd explicitly that the soft contributions are indeed small. We argued in particular 
that the contribution of collective modes is artihcially amplihed when the coupling is not 
too small by truncating this contribution at order g^. 

The approach that we have developed relies on a hierarchy of scales that emerges when 
the coupling is small. It is an assumption that the structure identihed at weak coupling 
survives when the coupling grows, e.g. as we lower the temperature. The fact that 
we hnd approximations where such an extrapolation works supports the validity of this 
assumption. But of course a complete check of the method can only be done by comparing 
with “exact” results, such as those provided by lattice techniques. Much can be learned 
also from detailed comparisons with the other approaches discussed earlier, namely hard 
thermal loop perturbation theory or dimensional reduction. There is hnally much to do 
also within the present theory itself, to push for instance non perturbative renormalization 
techniques (interesting progress in this direction has been reported recently |^; see also 

li)- 
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Further tests involve the calculations of different observables. An example, advocated 
by Gavai at this workshop, is that of quark susceptibilities. This is a calculation 
that we have recently taken up, with interesting results The calculation of quark 


susceptibilities represents a small incursion into the physics of hnite chemical potentials 
which is of course readily available in our approach [jT^ . 
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